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Abstract 



We make an explicit formulation for the proton decay rate in the minimal 

•T— I , 

^ \ renormalizable supersymmetric (SUSY) SO(IO) model. In this model, the 



Higgs fields consist of 10 and 126 SO(IO) representations in the Yukawa in- 



teractions with matter and of 10, 126, 126, and 210 representations in the 
Higgs potential. We present all the mass matrices for the Higgs fields con- 
tained in this minimal SUSY SO(IO) model. Finally, we discuss the threshold 
effects of these Higgs fields on the gauge couplings unification. 

PACS number(s): 12.10.-g, 12.10.Dm, 12.10.Kt 
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I. INTRODUCTION 



Proton decay would be a smoking gun signature for Grand Unified Theories (GUTs). 
Unfortunately, no such signal has been seen. In fact, very strong experimental limits have 
been set for this process, placing the minimal GUTs in a very precarious position. Su- 
per Kamiokande (SuperK) has set a lower limit on the proton lifetime in the channel p — > K^V 
as 

t(p ^ X+P) > 2.2 X 10^3 [years], (1.1) 

at the 90% confidence level [1]. This has already placed stringent constraints on SU(5). In 
fact, minimal renormalizable SUSY SU(5) model is almost absolutely excluded [2]. ^ Thus 
the realistic unified model builders must seriously consider the proton life time constraints. 

Now, SO (10) GUTs have been mainly discussed in connection with the neutrino oscil- 
lations since this part reveals the physics beyond the Standard Model. In this connection, 
SO(IO) GUTs have some advantages over SU(5) GUTs. One of them is that they incorpo- 
rate the right-handed neutrinos as the member of the 16 dimensional spinor representation 
together with the other standard model fermions, and provide the natural explanation of 
the smallness of the neutrino masses through the seesaw mechanism [4]. In this paper, we 
consider the minimal renormalizable SUSY SO (10) model. This model contains two Higgs 
fields 10 and 126 in the Yukawa interactions with matter [5] [6]. This is a minimal model 
in the sense that it contains only the renormalizable operators at the GUT scale and it has 
the minimal contents of the Higgs fields compatible with the low-energy experimental data. 
If we relax the renormalizability at the GUT scale, the different minimal SO (10) models are 
also possible to consider [7] [8]. In this paper, we restrict our arguments within the renor- 
malizable theory at the GUT scale, and use the name of minimal in this restricted sense. 



If we take some textures of the mass matrices for fermions and sfermions, we may get a safe 
region for the proton life time in a minimal renormalizable SUSY SU(5) model [3]. 
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As was shown in [5] [6], this theory is highly predictive. However, the recent data [9] [10] 
showed that one of our prediction about the neutrino mass square ratio , "^|°' ~ 0.19 [6] is 

atm 

out of 3(7 allowed region. So we need to make re data-fitting compatible with the up to date 
experimental data well. However, the development of GUTs and rich experimental data 
drive us to a new stage of precision calculation. That is, we must do the precise estimations 
to include the ambiguities, coming from the input data and the threshold corrections. The 
former comes from the fact that there are many input data and they have rather large error 
bars (see, the strange quark mass, for instance). We must scan these many possibilities 
systematically. This is tedious but rather technical and is on going in a separate form. The 
latter is more fundamental. It depends on the details of the superpotential of the Higgs sec- 
tor, whose effects are not confined in the low-energy data predictions. In order to investigate 
the proton decay rate and the gauge couplings unification in precise, we have to determine 
all the mass spectra of the Higgs fields in terms of the parameters presented in this model. 
This is a very complicated task itself and is the main motivation of this work. Even in the 
minimal model, there are so many free parameters. So, in practical analysis of the proton 
decay rate and also the gauge couplings unification, one has to reduce the number of free 
parameters. That means we should consider the smallest number of the Higgs contents. 
Thus, we introduce our Higgs system as the most simplest one, {lO© 126© 126©210}. 
The meaning of the introduction of these representations will be revealed in the next section. 
Since our results are the general one for the "minimal" renormahzable SO(IO) models, it 
can be applicable to any parameter regions. For instance, even if we fix the type of the 
Yukawa couplings in the matter sector and also the Higgs potential, the result is not unique. 
If we restrict the values of the parameters in the superpotential to some restricted region, 
we may get the two different types of the seesaw mechanism, type-I [5] [6] or type-II [11]. 
In this paper, we do not explain the way to save the model from the proton decay rate of 
SuperK, explicitly. Our main purpose of this paper is to produce all the mass spectra of 
the Higgs fields including all the Clebsch-Gordan (CG) coefficients and to propose a general 
formulation which is applicable to any parameter choices. In these applications, our theory 
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might be found to be insufficient. Even if it is the case, our theory is very useful for more 
elaborate theory. 

This paper is organized as follows. In Sec. II, we give the explicit form of the superpoten- 
tial in our model. In Sec.III, a very brief description of the symmetry breaking procedure 
and the decomposition of the original Higgs fields into the minimal supersymmetric standard 
model (MSSM) are given. In Sec. IV, using these techniques, we can get the mass matrices 
for a variety of fields, especially for the would-be Nambu-Goldstone (NG) modes. Then we 
can check that the appropriate NG modes do appear in the mass spectra. In Sec.V, we check 
the mass matrices for the electroweak Higgs doublets and consider the conditions for two 
Higgs doublets to remain light. In Sec. VI, we derive the formulae for the evaluation of the 
proton decay rate. In Sec. VII, we finally check the remaining mass matrices and the effects 
of the threshold corrections on the gauge couplings unification. In Appendices, we list up 
all the coefficients of dimension-five and -six operators, which are relevant to proton decay. 
The applications to a more elaborate model will be given in a separate publication. 

II. MINIMAL SO (10) GUT 

In this section, we explain the minimal renormalizable SUSY SO(IO) model. As men- 
tioned in the introduction, it contains two Higgs fields in the Yukawa interactions with matter 
[5] [6]. In the SO(IO) models, the left- and right-handed fermions in a given i-th generation 
are assigned to a single irreducible representation 16i = ^i. Since 16(8)16 = 10®120®126, 
the fermion masses are generated when the Higgs fields of 10, 120, and 126 dimensional 
representations develop nonvanishing vacuum expectation values (VEVs). The use of only 
one Higgs field, 10 in the Yukawa interactions with matter is obviously ruled out for the 
description of the realistic quark and lepton mass matrices. Furthermore, the use of 126 
dimensional Higgs field has desirable properties for providing masses of the right-handed 
Majorana neutrinos. Also it was found that 10 (= H) and 126 (= A) are suitable for the 
mass matrices since they satisfy the Georgi-Jarlskog relation. In order to preserve super- 
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symmetry, we must also include the Higgs field A of 126 dimensional representation. The 
Higgs field $ of 210 dimensional representation is introduced to break the SO(IO) gauge 
symmetry [12] and to make mix the Higgs doublets included in H and A [5]. Then the 
minimal Yukawa coupling becomes 

Wy = Yl^^ifiH^Ifj + nie^iA*.-, (2.1) 

and the minimal Higgs superpotential is [12] [13] [14] 

W = mi$2 + maAA + msH^ + Ai$^ + Aa^AA + X^'^AH + \4<^AH. (2.2) 

The interactions of 210, 126, 126 and 10 lead to some complexities in decomposing the GUT 
representations to the MSSM and in getting the low energy mass spectra. Particularly, the 
CG coefficients corresponding to the decompositions of SO{10) — > SU (3)c x SU (2) l x U{1)y 
have to be found. This problem was first attacked by Xiao-Gang He and one of the present 
authors (SM) [15] and further by Lee [13]. But they did not present the explicit form of 
mass matrices for a variety of Higgs fields and also did not perform a formulation of the 
proton life time analysis. In this paper we will complete that program in the frame of our 
minimal SO (10) model. 

III. SYMMETRY BREAKING 

In order to discuss the symmetry breaking pattern, here we briefly summarize our con- 
ventions for the SO(IO) indices. SO(IO) indices a = 1, 2, • • •, 9, are divided into two parts 
q; = 1, 2, 3, 4 for S0{4) ^ SU{2) x SU{2) and a = 5, 6, 7, 8, 9, for S0{6) ^ SU{4). For the 
SO{10) — > SU (3)c X SU {2)lxU (l)y decompositions it is very useful to define a "F diagonal 
basis" ; 1' = 1 + 2i, 2' = 1 - 2^, 3' = 3 + 4i, 4' = 3 - 4i, 5' = 5 + 6i, 6' = 5 - 6i, 7' = 7 + 8i, 
8' = 7 — 8i, 9' = 9 -I- Oi, 0' = 9 — Oi (up to normahzation factor, l/-\/2). Hereafter we use 
this Y diagonal basis and omit the dashes: The 10 dimensional irreducible representation, 
H is spanned by the states a = 1,2, • • -,9,0. The 210 dimensional irreducible represen- 
tation, $ and the 126 © 126 dimensional reducible representation A -|- A, are spanned by 
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the anti-symmetric tensors of the fourth rank {a(5^5) and the anti-symmetric tensors of 
the fifth rank [aP'-ySe), respectively. Here and below the bracket (■ ■ ■) represents the total 
anti-symmetrization of the indices within the bracket. 

The Higgs fields of the minimal 5'O(10) model contain five directions which are singlets 
under SU{?>)c x SU{2)l x U{1)y. The corresponding VEVs are defined by 

($) = ^0,4 (3.1) 

i=l 

{A)^vrv^, (A)^v^€^, (3.2) 
where 0j (i = 1, 2, 3) are included in 210, 

01 = -4| (1234) , (3.3) 



02 = ^ (5678 + 5690 + 7890) , (3.4) 

v72 



vr is in 126 



03 = -^^^ (1256 + 1278 + 1290 + 3456 + 3478 + 3490) , (3.5) 



-L= (13579), (3.6) 



and vr is in 126 



Vl20 



Notice that 



0i • (pj = = 1,2,3) , 

• = ^ • ^ = 0, 
v^-v^^l. (3.8) 

Due to the D-fiatness condition the absolute values of the VEVs, and vr are equal, 

l^^l = \vr\- (3.9) 



6 



Now we write down the VEV conditions which preserve supersymmetry, with respect to the 
directions 0i, 02, 4>3, and vr, respectively. 



(>2 



2mi0i + 3Ai-^ + A2^^ = 0, (3.10) 
6^ 10^6 ^ ^ 



^2 + 03 \ , X Vr-Vr 



2„,,,^3A.^-^J.A.^ = 0. (3.11) 
2m.^3 + 3A.(f^ + ^)+A,^^ = 0, (3.12) 

Here we consider only the solutions with \vr\ ^ 0. Eliminating vr ■ vr, (pi and 02 from Eqs. 
(3.10)-(3.13), one obtains a fourth-order equation in 03, 



(03 + ^) {8 0^ - 15 Mi0^ + 14 A^?03 - 3 A^? + (03 - Ml)' M2} = 0, (3. 14) 



where 



M, ^ 12 (^) , ^ 60 (^) . (3.16) 

Any solution of the cubic equation in 03 is accompanied by the solutions 

03 {Ml - 5 0i) 

01 = 



02 



(Ml - 03)2 ' 

1 (M2-2Mi03-0l) 
\/2 (Ml -03) 



_ 5 f\i\ 03(Mi- 3 03)(M2 + 0i) 

■ ^ 3 UJ (A^i-03)' ■ ^'-''^ 

The linear term gives the solution of the fourth-order equation (3.14) which is very sim- 
pie, 03 = -6 (^). It leads to 0i = -Vd (ff), 02 = -SV2 (^) and ^{vr-W) = 



(^) V ^ (m2 ) ~ (^) ■ '^^^^ solution preserves the SU (5) symmetry. Therefore, it is 
physically not interesting. The cubic term solutions lead to the true SU{?>)c x SU{2)l x 
U{1)y symmetry. 
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IV. WOULD-BE NG BOSONS 



In order to check the number of NG modes we write down the mass matrices for the 
Higgs(ino) fields which transmute the non-MSSM 5'O(10) gauge fields into very massive 
gauge fields. At first, we list the quantum numbers of the would-be NG modes under 
SU{3)c X SU(2)l X C/(l)y. 

(3,2,1)0(3,2,-1) 

(3,2,-i)©(3,2,i) 

(3,1,-1)0(3,1,1) 

(1,1,1)0(1,1,-1)], 

(1,1,0)]. 

Total number of the NG degrees of freedom is : 12 + 12 + 6 + 2 + 1 = 33. In the following 
subsections we give explicit expressions for the mass matrices and check that their deter- 
minants are zero. The mass matrices receive contributions from the F terms in the Higgs 
potential. The matrix elements of the mass matrices comprise the CG coefficients which 
appear as coefficients of the triple products of the SU{3)c x SU{2)l x U{\)y components 
of the Higgs superfields. For the calculation of the CG coefficients, one must first find the 
exphcit expressions for the SU{?>)c x SU{2)l x U{1)y components of the Higgs superfields. 
We will publish the complete tables of the CG coefficients of a more general Higgs sector in 
a separate publication and we will list only the mass matices in this paper. 

Note that the mass matrix for every irreducible representation under SU {3)c x SU (2)^ x 
U{1)y with Y ^ and the mass matrix for the corresponding complex conjugate represen- 
tation are equal up to transposition. Therefore, only one of the two accompanied mass 
matrices is listed. Of course, when enumerating the total degrees of freedom, one has to 
be careful to include all the mass eigenvalues (472 in total). The mass matrices define the 
mass part of the superpotential as a bilinear form of the fields and corresponding complex 
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conjugate fields. The basis for the mass matrix is defined as a row of the fields multiplying 
the mass matrix form the left. 



A. 



(3,2,1)0(3,2,-1) 



f (3 2—-) (3 2— -)1 

In the basis < $(6 2 2)'^ ' ^(10 2 2^) f (here and hereafter the lower indicies indicate SU {A)c x 
SU(2)l X SU{2)r and the upper SU{'i)c x SU{2)l x U{1)y in the case of double indices), 
the mass matrix is written as 



V 



2mi 



Ai<j!>3 
3\/2 



6 3V2 

9m, -I- _ Ai03 

^"^1 ^ 3V2 6 



(4.1) 



This determinant is indeed zero assuming the VEV conditions, Eqs. (3.10)-(3.13). 



B. 



(3,2,-i)e(3,2,i) 



(3 2 -) (32 i) (3 2 i) (3 2 i) 

In the basis i ^(6 2 2) ' ^(10 2^*2)' ^(15 2*^2)' ^(15 2''2) \^ mass matrix is written as 



2mi + ^ 
3V2 



V 



^2VR 
IOV3 





3t/2 



+ + 6 



5\/6 





X2VR 
10t/3 



1712 + 



^202 

30^2 





+ 



A2</'3 

20 







I A2(/'2 I A2 03 

"^2 + ^ + -60" / 



This determinant is also equal zero assuming the VEV conditions. 



(4.2) 



C. 



(3, 1,-1) ©(3, 1,1) 



(3 1 -) (31 -) (3 1 -) 

In the basis i $^.^5 ^^j^, $^-^g ^^3^, A^j^^ ^^3^ \, the mass matrix is written as 



2m, + 

Al</'3 

3\/2 

_ A2£R 



Al(/'3 

3a/2 
Ak/. 



2^1 + ^ + 11 



. A2£R 
IOV3 

5\/6 



A2f 



V lOv^ s'v^" + -T 30 y 

This determinant is also equal zero assuming the VEV conditions. 



A2<;^i 



+ 



A202 _j_ A2 03 



(4.3) 
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D. [(1,1,1)0(1,1,-1)] 



In the basis {^(15 13), ^[ini 3)}) mass matrix is written as 

t/2Ai02 X2VR ^ 

1712 + 



2mi + ^ + ^ 



10 



IOV2 / 



(4.4) 



10 ""^ ' 10V6 

This determinant is also equal zero assuming the VEV conditions. 



E. [(1,1,0)] 



In the basis | ^Ji^f,!)' ^(i^fi)' ^(^iJis)' ^(TO^^i) [' ^^^^ matrix is written as 



2mi 


\2VR 

X2VR 

\ 10a/6 



2mi + 



V2A1 



a/2Ai(/)3 

3 

_ A2£R 

10\/2 

_ 



2mi 



\/2Ai(^3 
3 

Ai(^i _|_ 
V6 + 

\2VR 
10 

_ A22fl 

10 



v^Ai 





A2ffi 


10^6 


lOv^ 


>^2VR 


>^2VR 


IOV2 


10\/2 


A217r 


A2ffl 


10 


10 


e4 








e4 



(4.5) 



Here e4 = 777,2 + A2 

( We + W2 + ^) nothing but Eq. (3.13) and Eq. (4.5) has one zero 

eigenvalue. 



V. ELECTROWEAK HIGGS DOUBLET 

In the standard picture of the electroweak symmetry breaking, we have the Higgs dou- 
blets which give masses to the matter. These masses should be less than or equal to the 
electroweak scale. Since we approximate the electroweak scale as zero, we must impose a 
constraint that the mass matrix should have one zero eigenvalue. 

We define 



0-10 _ t7-(l>2,i) 

-"u — -"(1,2,2) ' 



^(15,2,2) 



— ■^(15,2,2) 



- "^(10,2,2)- 



and 



(5.1) 
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M, 



doublet 



o-lO _ t7-(l'2-i) _ (1,2-i) _ (1,2,-i) _ {1,2,-i) , . 

-"d =-"(1,2,2) ' = ^(15,2,2) ' = ^(15,2,2) ' = *(10,2,2) • \P-^) 

In the basis A„, A^, the mass matrix is written as 

/ \ 

O™ A3 02 A3 03 \a4>2 A4(/)3 A4ffl 

A402 A403 I A202 A203 n A 

2v^ '"2 30 ^ ^ 

_ A302 _ A3 03 rv m _l_ ^2'/'2 _l_ A203 _ A22fl 

2V5 ^ "^2 + + 30 10 

As^lfl n _A2Uh I Al02 I Al03 

V5 ^ 10 zmi-h ^ + 2 / 

(5.3) 

The corresponding mass terms of the superpotential read 

Wm = A,, A,,$„) Mdoubiet Ad, 5^,$^)"". (5.4) 

The requirement of the existence of a zero mode leads to the following condition. 

det Mdoubiet = 0. (5.5) 

For instance, in case of A3 = 0, m2 + — = 0, we obtain a special solution to Eq. 
(5.5), while it keeps a desirable vacuum and it does not produce any additional massless 
fields. However, we proceed our arguments hereafter without using this special solution. 
We can diagonalize the mass matrix, Mdoubiet by a bi-unitary transformation. 

U* Mdoubiet = diag(0. Ml, M2, M3). (5.6) 

Then the mass eigenstates are written as 

hi, hi h^) = ,A„,A„,$J U^, 
{h,, hi hi h^) = {hH Aa, A,, V^. (5.7) 

The representations 45 and/or 54, and higher dimensional operators, are not included in 
our minimal model. Therefore, we must set the " Doublet-Triplet sphttings" by hand as Eq. 
(5.5). In the case of A3 = A4, Eq. (5.3) becomes symmetric, and Hu and have the same 
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coefficients in Eq. (5.7). Tfiis can not be accepted since it leads to tlie formal singularity in 
the low-energy Yukawa couplings (A matrix in Eq. (6.14)). Namely, it leads to the equality 
Yu — Yd, and therefore only the ratio of Fio and Y^e can be determined from Eq. (6.14). So 
we set A3 7^ A4 hereafter. 

By making the inverse transformation of Eq. (5.7), the following expressions are obtained, 

Hl^ ^auH^ + ---, Hf = aaHa+---, A„ = + • • • , 'K^^ (3dHa + ■ ■ ■ , (5.8) 

where "H " represent the heay^ Higgs fields, h^^ [i — 1,2,3) which are integrated out 

when considering the low-energy effective superpotential. 
Precisely, we can read off from Eq. (5.7) as 

Oin = (C/*)ll, = (C/*)l2, O^d = (l^*)ll, = (^*)l3. (5.9) 

Using the two pairs of the Higgs doublets, and A^^^, the Yukawa couphngs of Eq. (2.1) 
are rewritten as 

Wy = u1 {Yil Hi' + Yii, A„) + d1 [Yii Hf + YH, A,) 
+ {yll K - ^ylk A.) + 4 {Yll Hf - 3Yli, A,) £, 
+ (yiLvn) (5.10) 

By using Eq. (5.8), we obtain the low-energy effective superpotential which is described by 
only the light Higgs doublets Hy_ and if^, 

VFeff = u\ (a^Yll + Z^^nie) + d\ (aaYH + (5^X11^ qj 
+ {a^Yii - SPuYiL) Huij + el - Zf^aYH,) Hdij 

+ liesH^Hd. (5.11) 

Here we have assumed that some mechanism, like the Giudice-Masiero mechanism [16] in 
supergravity, may produce the effective /i term, /Xeg for the light Higgs doublets. 



12 



VI. PROTON DECAY 



After the symmetry breaking from SO(IO) to SU{3)c x SU{2)l x the generic 

Yukawa interactions between the matter fields and the color triplet Higgs fields are given by 



Wy = Yl^ (qdj + Wrd^) + Yli, {q^i, + u^.d'^) 



(6.1) 



Here we have defined 



TT r7-(3>l>3) TT T7-(3>1> 3) A 

-"r = -"(6,1,1) ' = -"(6,1,1) ' 



A (3,1,5) A _ A (3,1,-5) A ■ _ A V",-,-3 

^(6,1,1) ' '^T = ^(6,1,1) ^ ^-r = ^ 



(3,1,-i) 
(10,1,3) 



(6.2) 



For later use we define 



A(3,l,i) A 

^(6,1,1) ' 



(3,1,-5) A/__ a(3,1,3) ^_ 
' — ^(T0,l,3) ' 



AV"'"' 3 

^(6,1,1 



^(3,1,-i) 
^(15,1,3) ■ 



(6.3) 



In the basis A^, A^, A^j, the mass matrix reads 



-^triplet = 



2m3 



A4(/'l 



A3(/'2 

30 



VTo 



V30 



A3^>il 



15 



A4</'1 A4(/)2 



_ A2£R 

10a/3 

A2</'3 

15^2 



A3 01 I Xz(t>2 

7; ~r 



a/IO ^ V30 



m2 





A4t^i?, 

lOv^ 


m44 

_ A2^ 



where ^44 = 2mi + ^ + + ^Ai^ and mgs 



™ I A201 I X2<t>2 



The corresponding mass terms of the superpotential read 



\/2A4</>3 

a/15 

A2(/'3 

15^2 



5t/6 



(6.4) 



= (i/^, A^, A^, A^) MtHpiet (i^r, At, At, $t, A^)"^ . 



(6.5) 



Here we integrate out the color triplet Higgs fields. At and $t, which do not appear in the 
Yukawa interaction with matter, Eq. (2.1), 



13 









1 







2mi 



V6 



l4>2 I 2Al(/'3 A / ^3(^1 _|_ A302 ^ rj 

;V2"^ 3 J ' \ VTo ^ VsoJ 



X 

3V2 



1712 ^ 



^2 5^6 



(6.6) 



where P = ■ (2mi + ^ + ^ + 2^). 

Putting this into the original mass terms of the superpotential Eq. (6.5), we can obtain 
the following mass terms for the color triplet Higgs fields, 

= A^) M^i^, (Ht, At, A^)^ . (6.7) 

{rriij} are given as 



Here the explicit forms of the elements of this mass matrix, M^^p^^ 
follows. 



mil = 2m3 - — 



X3VR XiVR 

H 7E- ■ • 



A401 ^ A402 



30 



\rn + '^^'^^ + '^^'^^ _^ 2Ai03' 
, a/6 3\/2 3 ^ 



^301 _^ ^302 



10 



mi2 = 



A401 



A402 



10 



30 2^ 10^3 



m2 



X4VR 



mi3 
m2i 

m22 = 1712- 

m23 



V^X4(f)3 , 1 A2VR X4VR 



15 V 

A301 A302 



10 V30 

1 A2fiJ 



]_ X3VR X2V^ 

V V5 '"^^ ' 10^3' 

A2^^ 



m2 



m32 
m33 



V lOx/3 

A203 _ J_ A2t'R 

15^/2 2^ IOV3 

\f2X-i(p^ , 1 Agt'^ 



10^3' 

A2^^ 



m2 



15 V5 

A203 1 A2fR 



• m2 ■ 



15^2 10^3 

A201 A202 
777.2 + ^ + 



r?72 • 



10a/3' 

A2^ 

5^6' 

A2^^ 

5v/6' 
I A2t;ij 



10^6 ' 30V2 V 5^6 '^^ 5^6' 
Moreover, integrating out the color triplet Higgs field A^, we obtain the effective Yukawa 
interactions between the matter fields and the color triplet Higgs fields as 

Wy = % [qdj + u^d'^) + A^ (qdj + u^df) 



X2VR 



(6.8) 



+ ^10 Ht -QiQj 
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m33 

+ ^126 TiQiQj 



126 



2 

+ (l - ^) Yli, At {ute^ + d^u^) . (6.9) 
Then the effective mass terms for the remaining color triplet Higgs fields are written as 

^Hj.(aHT + bAT) 



+ AYicHr + dAr) 

= [h^, Aj) Mt 



(ht\ 



(6.10) 



where a, 6, c, d are defined by 



nils 

a = mil rusi, b = rriu ^32, 

m23 , m23 
c = m2i ■ msi, d = 17122 ■ m32. 



(6.11) 

m33 m33 
Combining the Eqs. (6.9) and (6.10) leads to the effective dimension-five interactions after 
integrating out the remaining color triplet Higgs fields [17], 



(6.12) 



inducing the dangerous proton decay. Here, Cl and Cr are given by the Yukawa coupling 
matrices at the GUT scale, Mq 



Cf\MG) = {Yil Yii,) 



1 ^ 



10 



vkl 

y -'^ 126 y 



Cf\MG)^[Y, 



10 



mi3 
m33 



-'126? 



^32 \ 
^33/ 



^126 I ^'^^r ^ 



10 
vkl 

y -"^ 126 y 



(6.13) 



Note that 



Y, 



10 



yll26y 



V" (v\ 



Y, 



— A-l 



= A 



(y\ 



15 



(6.14) 



Thus we have 



( yki \ 

u 
Vkl 



(6.15) 



We make use of this expressions in order to evaluate the renormalization group effects on the 
Wilson coefficients Cf^'^'' and C]|*^^ Without loss of generality, we can use the basis where 
Yu is real and diagonal, 

y„ = — ^—T diag(m„, rric, m*), (6.16) 
vsm p 

with V ~ 174.1 [GeV]. Since is a symmetric matrix, it can be described as 

Yd = ^ ^CKM diag(md, m^, m^) T^ckm> (6-17) 

I'COSp 

by using a unitary matrix 

FcKM = e*"i e*"^^^ e*"3^« Vckm e'^''^' e'^'^' , (6. 18) 

where A3, Ag are the Gell-Mann matrices and Vckm is the Cabibbo-Kobayashi-Maskawa 
(CKM) mixing matrix [18]. ^ 

The complete anti-symmetry in the color indices requires that the dimension-five operator 
Eq. (6.12) possesses the flavor non-diagonal indices [19]. As a consequence, the dominant 
decay mode is p — > K'^P. This fact implies that the chargino dressing diagrams dominate 
over the gluino and the neutralino dressing diagrams [20] . 

In the components form, the dimension-flve operators at the SUSY breaking scale, Msusy 
are written as 

£,^(udue)XYij T' , /~i(uude)XYij^ — j 

5 = Q 'uxdYULieLj + Cl ^uxUrdueLj 

, y^(udue)XYij , ^(uude)XYi'j ^ r 

+ C)j uxdyURieRj + C)j ■ -uxUvdRiCRj 

I /^(uddv)XYij T' 7 . ^(dduv)XYij^ 'T' j" 

+ Cl UxdydLil^Lj + Cl -dxdYULiVLj 



^In Ref. [6], we set these phases {i = 1, 2, 3), (3i {i = 2, 3) to zero or tt. 
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+ Q UxeyULidLj + Cl -dxeyULiULj 



y{deuu)XYij 1 

I ^(ueud)XYij — j . j^(deuu)XYij 1 T' — 

+ C)j uxeyURidRj + -dxeyURiURj 



^{deuu)XYij 



^{duud)XYij, ~ , ^{ui'dd)XYiji-~~ , , 



(6.19) 



The coefficients are obtained from tfie coefficients of tfie original dimension-five operators 
including their renormalizion from Mq to Msusy- Their explicit forms are found in Appendix 
A. After the sparticles dressing, we obtain the following type of dimension-six operators 
causing nucleon decays, 

1 



+ C'j^r^^'\uLdLi){uReRj) + C^RR''^'\uRdRi){uReRj) 

+ C^ll'^''^^^'^ {'^LdLi){dLjVLk) + Cr^'^"^'^'' {uRdRi){dLjl^Lk) 

+ 4l 



167r2 

~i{udue)ij I 



^i'^'^^^)^^^hdRdR,){uLl^Lk) 



(6.20) 



These operators should be renormalized from Msusy to Mz and further to the hadronization 
scale (//had) ~1 [GeV]. Then the effective four- Fermi Lagrangian is converted to a hadronic 
Lagrangian by using the chiral Lagrangian method [21] [22]. Details are given in Appendices 
B and C. 

For the decay mode p — > K^V^, the partial decay rate is given by the formula 



Tip ^ K+Vi) 



327r 



%1 -^lAip^K^u, 
ml ) U 



(6.21) 



Here nip — 0.938 [GeV] is the proton mass, mK+ — 0.493 [GeV] is the kaon mass and 
= 0.131 [GeV] is the pion decay constant. 
The amplitude A{p K^Vi) for p — > K^Vi reads [23] 

2mr 



A{p K+ 



+ 



n/^{uddy)21i . ^{uddy)21i 
n^{uddv)12i . ^{uddu)12i 



3m£ 
1 + 



mp 

2)1713 



{3F + D] 



, ^(ddui^)l2i 

+ (^Crl 



nir. 



3m 



B 



(3F - D) 



(6.22) 
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Here niB — 1.150 [GeV] is an averaged baryon mass, F — 0.44, D — 0.81 are the parameters 
in terms of which the octet-baryon axial- vector form factors are expressed and a, (5 are the 
hadron matrix elements which are defined by [24] 

Q;tiL(k) = (0|(iijtiRtiL|p(k)), 

(3ul{\^) = {Q\dLULUR\p{\^)). (6.23) 

Wi:(k) denote the left-handed components of the proton wave function. It is known that 

I a I = 1/3 1, and {3 is in the range [24] 

0.003 [GeV^] < /3 < 0.03 [GeV^]. (6.24) 
From the recent lattice calculations, one group reported that [25] 

q; = -(0.015 ±0.001) [GeV^], 

/3 = 0.014 ± 0.001 [GeV^] . (6.25) 
But the other group reported the smaller values [26] 

a = -(0.006 ±0.001) [GeV^], 

(3 = 0.007 ± 0.001 [GeV\ (6.26) 



VII. GAUGE COUPLINGS UNIFICATION 

In general, the gauge couplings unification imposes constraints on the mass spectrum of 
many varieties of Higgs fields [27] . Our strategy is a generic one that all of the dimensionless 
coefficients should remain of order one to preserve the perturbative limit and put all the 
VEVs at the GUT scale in order to realize the simple gauge couplings unification picture. 
For the numerical evaluation, we use the one-loop renormalization group equations (RGEs) 
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in the DR scheme [28] , 



3 4 



1 1 , C2{Gi) , 1 



+ 



MgJ y ' VaC"^""'^ 



(7.1) 



where C2 is the quadratic Casimir operator; C2 {SU{3)) = 3, C2 (-SC/(2)) = 2, C2 {U{1)) = 0, 
and C denotes the Higgs fields which have the corresponding gauge quantum numbers. 
is it's mass matrix and "dct' " means that the determinant should be taken excluding the 
zero modes, bi and 6^ are the /? function coefficients ; 63 = —3, b2 — 1, bi — and b^ are 
given in Tables I and II. For ai (Mz) l^s, we use the following values. 

«3 (Mz) Ims = (Mz) , (7.2) 



«2 (Mz) Ims = « (Mz) / sin^ (Mz) , (7.3) 
5 
3 



ai (Mz) Ims = ^ « (^2) / (l - sin' (Mz)) , (7.4) 



with [29] 



(Mz) = 0.1172, a (Mz) = 1/128.92, sin' {Mz) = 0.23113. (7.5) 

Excluding the fields which mix with the would-be NG fields and the fields with SU{3)c x 
SU{2)l X U{1)y quantum numbers, [(1,2, |) +h.c.] ("Higgs doublet") and [(3, 1, |) +h.c.] 
(color triplet Higgs fields), the massive fields are given as follows. 



For 126 and 126 representation fields, their quantum numbers, the masses and their P 
function coefficients are given in Table I. 

For 210 representation field, their quantum numbers, the masses and their (3 function 
coefficients are given in Table II. 



DR uses dimensional regularization through dimensional reduction with modified minimal 
subtraction. 

^Here we assume, for simplicity, all the mass eigenvalues of the Higgs fields are smaller than Mq 
and all the masses of the gauge fields lie around Mq- In the other cases, the formula becomes quite 
complicated. 
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TABLES 



TABLE I. The mass matrices and the /? function coefficients for 126 and 126. 



quantum numbers 


mass matrices, or mass eigenvalues 










(8,2,^) 


+ h.c. 




^ rn-^ '^2</'2 A2c?i3 n ^ 
"^2 30V2 60 U 

n mr, '^202 1 A203 




12 


8 


24 
5 




(6,3,i) 


+ h.c. 


^^"2 lov^ 30v^ 


15 


24 


12 
5 




(6,1,1) 


+ h.c. 


1 \24>1 ^2<t>2 ^2<t>S 

'^"''^ ^ IOV6 30v^ 30 


5 





5 




(6,1,1) 


+ h.c. 


^ 1 X2<i>2 1 A2</'3 

10V6 30a/2 ' 30 


5 





1 a 
5 




(6,l,i) 


+ h.c. 


I A2<jil A2<^2 

"''^ ^ 10^6 30^2 


5 





4 
5 




(3,3,i) 


+ h.c. 


A201 1 A2<?!-2 

"^2 10V6 ' 30V2 


3 


12 


6 
5 




(3,2,1) 


+ h.c. 




( mr, 1 ^2</'2 A2</>3 A ^ 

"^2 30V2 60 ^ 

n ™ 1 A2</>2 A2</>3 

^ " "^2 + 20 ) 




2 


3 


49 
5 


(3,1,1) 


+ h.c. 


™_ 1 A2<pl 1 A2<^2 A2<^3 

"''^ ^ IOV6 ' 30V2 30 


1 





32 

5 


(l,3,l)+h.c. 


A2</>1 I A2<?i2 





4 


18 
5 


(l,l,2)+h.c. 


1 ^2(I>1 1 A2<^2 A2</'3 

^^"2 ^ 10^6 ' 10a/2 10 








24 
5 



20 



TABLE II. The mass matrices and the (3 function coefficients for 210. 



quantum numbers 


mass matrices, or mass eigenvalues 


4 


4 




(8,3,0) 


'f^ \t 


9 


16 





(8,l,l)+h.c. 




6 





48 
5 


(8,1,0) 




I'm . '^l'^2 Ai03 \ 
Al03 0™, 1 Ai0i Ai(^2 

V 3V2 ^""^ ' V6 3V2 / 










(6,2,1) +h.c. 


0^, Al02 Al(/)3 


10 


6 


10 


(6,2,i) +h.c. 




10 


6 


2 
5 


(3,3,1) +h.c. 




3 


12 


24 
5 


(3,1,§) +h.c. 


^^^^ VE ^ 3V2 3 


1 





10 


(1,3,0) 


^"^1 1 3 





2 





(1,2,1) +h.c. 


2mi + ^1/3 





1 


27 
5 
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Putting these values into Eq. (7.1), the unification condition produces two individual 
equations, 

as {Mg) = a2 {Mg) , (7.6) 

and 

as (Mg) = a, {Mg) . (7.7) 

Setting all VEVs at the GUT scale, 0i ~ 02 ~ 03 ~ \vr\ Mg, and the remaining 
dimensionless coefficients of order one, we can search whether Eqs. (7.6) and (7.7) have a 
solution for Mq below the Planck scale, Mq < Mpianck- If such a solution exists, it would 
limit the parameters in the superpotential Eq. (2.2) to some restricted region. 

VIII. CONCLUSION 

We find the general formulation for the proton decay rate in the minimal renormalizable 
SUSY SO(IO) models. Using this generic formulation one can find whether the minimal 
SUSY SO(IO) grand unified theory has been excluded. 

Recently, using their Yukawa couplings (Eqs. (8) and (9) in Ref. [30]), Goh-Mohapatra- 
Nasri-Ng obtained the allowed region of {x,y,z) which correspond to §) in our 

notation. However, they did not discuss the concrete form of the superpotential and, there- 
fore, compatibilities of their superpotential with the other constraints are not clear in their 
paper. Also, as we have mentioned above, there appears a non zero x value even without the 
54 dimensional Higgs field. Further, besides the color triplet Higgs fields, there is a much 
richer Higgs particle contents. These additional Higgs fields may cause a pathology of the 
gauge couplings unification. This paper presents a relationship among these comprehensive 
but tightly connected problems. 
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NOTE ADDED 

In the recent paper [arXiv:hep-ph/0402122] and the revised version of the paper 
[arXiv:hep-ph/0204097], the authors claimed that they disagree with our results for mass ma- 
trices. We point out that our results satisfy all possible consistency checks. Namely, for arbi- 
trary couplings mi^2,3 and Ai^2,3,4 in Eq. (2.2) there are solutions characterlized by = 0, 
particularly SU{5) x f/(l), SU{5) x U{1) flipped, SU{3)c x SU{2)l x SU{2)r x U{1)b-l and 
SU{3)c X SU{2)l X U{1)r x U{1)b-l symmetry preserving vacua, (see, Eqs. (3.10)-(3.12) 
with = 0). In all the above four symmetry breakings and in all the mass matrices, 
the field $ decouples from the other set of fields, H, A and A. Moreover, all our mass 
eigenvalues for the multiplets coming from the field coincide with the corresponding re- 
sults in [15]. Furthermore, for the SU{5) symmetry breaking case, \vr\ ^ (see, the end 
of Sec. Ill), all our mass eigenvalues (472 in total) and the corresponding multiplets under 
Sl]{2i)c X SU{2)l X U{1)y are grouped according to the SU{5) irreducible representations 
with the correct would-be NG fields. The details will be published in a separate publication. 

APPENDIX 

A. Dimension-five operators 

In this appendix, we list the explicit form of the various interaction coefficients. 
We use the following notations for the mixing matrices which diagonalize the squark, 
slepton mass-squared matrices and chargino, neutralino mass matrices. Squark, slepton 
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mass-squared matrix M~, chargino and neutralino mass matrices Mq and are diagonal- 
ized by the unitary matrices Uj, O^, Or and On, respectively. 

UjMfui=di^g{ml,....,ml), 

Or Mc 0\ = diag(m~- , m~- ) , 

O^MjvOjv = diag(m~o,m~o,m~o,m~o). (^.1) 
For the dimension-five operators, we have the following expressions ^ 

^iui.e)XYii ^ ^^.m _ Cf^'){m)x,,+,{W~)y,^,, {AA) 

^iuMe)XYr, ^ ^^.^ _ C-^^'){U^)x,k+s{U^)Y,W, (A.5) 

^iudd.)XY^J ^ ^^^nkl _ ^ Xk{U*j)Yl{VcKMU{UMNs) jn, {A.6) 

^(A,uu)XYi, ^ cf^\ui)xAU~)Yu (^.9) 

^(.e-.d)xy., ^ ^^.^m _ ct'''){lJi)xMm)y.i^^^ (^-10) 

^iMuu)XYr, ^ ^^.^m _ Ct'^){U*~)x,k+s{Uj)y,i^„ {A.ll) 

^(ii>ud)XYii ^ ^^kUm _ Cf''){U*~)xkmYl{VcKM)jm, {A.12) 

^(uudd)XY^, ^ ^^nlkm _ Cf'^){U^)xkmYl{VcKM)im{VcKM)jn. {A.13) 

In Eqs. (A.6) and (A. 7), it should be noticed that the neutrinos in the final states should 
be rotated from the fiavor eigenstates to the mass eigenstates by using the Maki-Nakagawa- 
Sakata (MNS) mixing matrix [31], C/mns- 

^We use a notation for an anti-symmetric tensor, A^'^^^^^ = A^^^^ — A^^'^K 
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B. Sparticles interactions 



We use the following notations for the quark-gluino-squark, quark(lepton)-chargino- 
squark(slepton) and quark(lepton)-neutralino-squark(slepton) interactions, 



quark-gluino-squark interactions 



gux — iV2d^ 



Gfi'^Pr. + Gfi'^Pj gdx + h.c. 



quark(lepton)-chargino-squark(slepton) interactions 



[B.l) 



XAdx + d1 



^L{d) p ^-R(d) p 



XaUx 



r^Lie) p ^R(e) p 



Xa^x + h.c. 



(5.2) 



quark(lepton)-neutralino-squark(slepton) interactions 



^^iAX^L + J-^iAX ^R 



x\ux + d. 



N^/JPl + nS^Pr 



+^^KSPR&x + e^ 



NfiiPL + Nfi^^Pn 



x\ex + h.c. 



(5.3) 



Explicitly, we have the following expressions 



G^x^ = g3{Ul)x,i+3, 



Gfx"^ = g^iUDxi, 



G 



L{d) 
iX 



9sm 



d)X,i+3, 



Gfx^ = g3{Ui)xk{VcKM)ik, 



-iL{u) ^ ruu^ 
^iAX ^^My^sin(5 



{0*^)A2mxi. 



(5.4) 

(5.5) 
(5.6) 
(5.7) 
(5.8) 
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1 



(5.9) 
(S.IO) 

(Rll) 

(5.12) 
(5.13) 

(5.14) 

(5.15) 

(5.16) 
(5.17) 



-^iAx = [(<^^)a2 - tanew{Olf)Ai] {Ut)x,k{UMm)ik, 



[B.ll , 
(5.19) 



iV^iJ ^ | ^J'J;3^ (Q^)a3(^I)x. + ltan9w{0*r,)Ai{UI)x,+3^ , (i^.20) 



2 
3 



(5.21) 



'"^'^ ^/2 1Mwcos/? 

These expressions are found in [32], but only for the quark sector. So here we write them 
explicitly. 
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C. Dimension-six operators 



For the dimension-six operator, we devide the coefficients into three parts according to 
the dressed sparticles, 



^{udue)ij ^{udue)ij 



{g) + Ctt^'^x') 



{C.l) 



etc. Then we have the following expressions. These expressions have the same forms as [23]. 
However, ours are different from them in the neutrino sector as was mentioned in the end 
of Appendix. A. 



^{udue)ij ^ 



3m~ 



{udue)XYlj ^R{u) ^R{d) jp 



2 2 
9 9 



2 ' 2 
Ux dy , 



{C.2) 



m 



^{udue)XYlj ^R{u) ^R{d) rp 



2 2 

m~+ m~-, 

'X-A 

2 ' 2 
Ux dy . 



+ c 



{dvud)XY\i^R{d) ^R(,e) 



lAX'^jAY 



171-+ m~ 

P I ^A ^A 



1 1 2 

dx . 



(C.3) 



^{udue)ij / ~0\ 



Mudue)XYlj ^jR(u) ^rR{d) ^ | X^i 

ux dy , 



. ^(ueud)XYli , 7-iJ((i) , 7-JJ(e) j-, 
+ <^L ^^lAX^^jAY^ 



2 2 

Xa Xa 

2 ' 2 



^{uduejij ,~\ _ 4 1 ^{udue)XYlj ^L{u) QL{d) p 



c 



C 



{udue)ij f ~o 



3 m- 



m 



^{udue)XYlj ^L{u) ^L(d) p 



Xa 



^{uduejXYlj jyj-L{u) jij-L{d) j-, 



2 

mA 

9 


9 


2 ' 2 1 

m~ m~ 

Ux dy / 


I Xa 


m~+ 
Xa 


\ m~ 

\ Ux 


' m~ 




Xa 


\ m~ 

\ ux 


' 2 

m~ 

dy 



(C.4) 
(C.5) 



. ^{ueud)XYli i^jR{d) i.T--R(e) t-, 



lAX-^^jAY' 



rriK mK 

Xa Xa 



2 ' 2 
Ux ey ^ 



c 



4 1 

3 m~ 



C 



{udue)XYlj^R{u) ^R{d) 



2 2 



R 



^nyu) ^ixyu) rp 
'^IX '^iY ^ 



2 ' 2 

m~ m~ 

ux dy , 



{C.l) 

(as) 
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m-. 



-C 



{udue)XYlj ^R(u) ^R{d) 



R 



1 1 

m~+ m~- 



(dvud)XYli^L(d) ^L{e) j-, 



1 1 

m~+ m~ , 

Xa Xa 

2 ' 2 

m~ m~ 

dx 



C 



{udue)ij / ~o 



LR 



^{udue)XYlj AT-R(«) Ar-R(d) 



lAX-^^iAY 



{ueud)XYli ^rL{d) i^Me) 



^^lAX^^jAY^ 



Xa 



2 ' 2 
dx 



c 



{udue)ij 
RR 



(9) 



4 1 
3m~' 



'R ^IX ^iY 



c 



{udue)ij I ~± 



RR 



(r) = 



-c 



{udue)XYlj ^L{u) ^L{d) 



^lAY^iAX 




2 ' 2 

m~ m~ 

Ux dy . 



c 



{udue)ij i~Q 
RR 



C 



{udue)XYlj t.tL{u) n,TL{d) rp 
^^lAX^^iAY^ 



Xa x "a 

2 ' 2 



+ C 



(ueud)XYli ^ A7-i?(e) 



R 



^A ^A 



ux 



2 



ey . 



^{uddu)ijk 



4 1 
3 m~ 



C 



(uddv)XY jk ^R{u) ^R[d) 



2 2 
9 9 



'IX 



lY 



1 ' 2 

m~ m~ 



, ^{dduiy)XYlk^R{d)^R{d) 



2 2 
9 9 



2 ' 2 

m~ m- 



dx 



dy. 



C 



{uddv)ijk , ~± 



LL 



'^x\ 



^(uddv)XY jk ^R{u) ^R{d) j-, 
^lAY^iAX^ 





2 






2 


2 




m~ 


Ux 


ay 



+ c 



(ue«(i)Xyii^il(u)^il(e) ■ 



2 2 
2 ' 2 

m~ m~ 
dx , 



^(u(ic^^/)^jfc^~0^^ 



^{uddu)XYjk ArR(u) jKTR{d) jp 
^^lAX^^iAY^ 



Xa r\ 



2 ' 2 



+ C 



idduu)XYlk ArR(d) ArR(e) 



jAX-^^iAY 



2 2 

Xa x\ 



+ C 



(di^ud)Xyii A ril(d) A7-ii(e) 



AT-n-W AT- 

-'^jAX-'^ifcAy 
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c 



{uddv)ijk I 
RL \ 



4 1 
3 m~ 



C 



{uddv)XY jk ^L{u) ^L{d) j-, 



1 2 



C 



{uddv)ijk I ~± 
RL 



(r) = 



m 



1 ' 2 

m~ m~ 

UX dY 

,2 



-c 



(uddu)XY jk^L{u) ^L{d) 



+ C 



{ueud)XYli^R{u) ^R(e) 



R 



^jAX^kAY^ 



'^lAY'^iAX-^ 



2 2 
2 ' 2 



Xa Xa 



2 ' 2 

m~ m~ 



dx 



^{uddu)ijk ( ~n\ _ 1 ^{uddv)XYjk t^tL{u) j^rLid) rp ( ^Xa 
^RL yX. ) — ~ ^^lAX^^iAY^ \Z:J2~-' 



C 



{dduv)ijk 
RL 



\9) 



4 1 ^{uddv)XY\k ^L(d) ^L(d) j-, 
'^r (Jix ^jY ^ 



Xa 



2 ' 2 

m~ m~ 

Ux dy , 

2 2 

9 a 



3 m~ 



2 ' 2 
dx dy . 



c 



(dduv)ijk / ~± 
RL 



c 



{dduv)ijk f ~0 
RL 



1(1 1 
X°A X°A 
ml ' ml 
dx dy . 



Here we have defined a loop function, 



(ai6) 

(C.17) 



(CIS) 

(ai9) 

{C.20) 
(C-21) 

(a22) 



F(x,y) 



X — y I 1 — X 



logx 



logy . 



(a23) 
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